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D. FRETTLOH 

Abstract. The concept of star-duality is described for self-similar cut-and-project tilings 
in arbitrary dimensions. This generalises Thurston's concept of a Galois-dual tiling. The 
dual tilings of the Penrose tilings as well as the Ammann-Beenker tilings are calculated. 
Conditions for a tiling to be self-dual are obtained. 



1. Introduction 

Nonperiodic discrete structures like infinite words, tilings, or point sets, are a rich source of 
exciting objects studied within many fields of mathematics, including automata theory, com- 
binatorics, discrete geometry, dynamical systems, group theory, mathematical physics and 
number theory. Many interesting examples of such structures can be produced by substitu- 
tions as well as by cut-and-project schemes. For such structures there are several concepts of 
'dual' tilings (words, point sets) [17], [21], [24] • The present article describes how Thurston's 
concept of Galois-duality [23] generalises to arbitrary cut-and-project sets. This generalisation 
does not depend on Galois conjugacy but on the star map (see [15], or Definition 11.51 below) 
wherefore the term star-dual (resp. *-dual) is used, rather than Galois dual. The benefits 
of star-duality are that it gives a purely algebraic approach to certain geometric problems; 
that it covers a pretty wide range of examples; and that it may help to friend the certain 
concepts of duality. Furthermore, and hopefully, it can be utilised to attack open problems 
in the theory of non-periodic structures. 

The first section of the present paper collects necessary facts and definitions about substi- 
tutions and cut-and-project tilings, resp. model sets. Section [2] contains the definition of 
★-duality, together with an example. The calculation of the *-dual tilings of the famous Pen- 
rose tilings and Ammann-Beenker tilings [TT] is carried out in Section [3l It turns out that 
the star-dual tilings of the Penrose tilings are the well-known Tubingen triangle tilings (and 
vice versa), which coincides with a related result in [2]. In contrast, the star-dual tilings of 
the Ammann-Beenker tilings are very similar to the Ammann-Beenker tilings themselves. In 
fact, they are mutually locally derivable [4] with the Ammann-Beenker tilings. This moti- 
vates the question whether there are self-dual tilings w.r.t. star-duality. Therefore, Section [4] 
is dedicated to self-duality in arbitrary dimensions. A necessary, and a sufficient condition 
for a tiling to be self-dual w.r.t. star-duality are obtained. 

The concept of star-duality is not new, see for instance [23], [10], [24J. The relevant results 
of this paper are the following: First, Theorem 11.31 gives a useful technical result which is 
exploited in the next sections, and will be exploited further in future work. Second, the dual 
tilings of the Penrose tilings and the Ammann-Beenker tilings are obtained, where the latter 
is done for the first time here, up to knowledge of the author, where the former one is obtained 
in [2] with much more effort in a slightly different framework (but not in [ID], as the title 
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Figure 1 . The substitution for the Penrose tiling in the version with triangles 
S and L (left), the first iterates of the substitution on the obtuse triangle L 
(right). The substitution factor in this example is r = v/ ^ +1 . 

might suggest). Last, Lemma 14.21 and Corollary 14.31 are first steps towards a classification of 
self-dual tilings, in one dimension as well as in arbitrary dimension. 

Remark: For convenience, we will exclude periodic structures in the following. Each occur- 
ring biinfinite sequence, tiling etc. is assumed to be nonperiodic. The set of positive integers 
is denoted by N. The closure (interior) of a set A is denoted by cl(A) (int(A)). 

1.1. Substitution sequences. Let A be an alphabet, i.e., a finite set of letters; let A* be 
the set of all finite words over the alphabet A, i.e., all words which are concatenations of 
letters of A. Usually, the empty word e (consisting of no letters) is required to be an element 
of A* . A substitution is a map ( : A — ► A* \ {e}. We will call such a substitution also 
symbolic substitution if we want to emphasise the difference to a tile-substitution. By setting 
C(ai • • • a m ) = C(ai) • • • C( a m), C ac ts also as a map from A* to .A*, from to A^, and from 
A z to A z . The family of all biinfinite words, where each finite subword is contained in some 
a k (a) (a £ A), is denoted X a . There is a lot of literature devoted to symbolic substitutions, 
see for instance |16] and references therein for more details. 

1.2. Substitution tilings. A tile is a nonempty compact subset of M d which is the closure of 
its interior. A tiling of M. d is a set {T^.^ of tiles which is a covering ofR d (i.e., M. d = U jeN Tj) 
that is non-overlapping (i.e., int(3l-) Dint^) = for j ^ i). A simple way to generate tilings 

- periodic ones as well as nonperiodic ones — is via a tile-substitution. Roughly speaking, 
a tile-substitution is given by a finite set {T ± , . . . ,T m } of tiles, the prototiles, an expanding 
map Q, and a rule how to dissect each QT- into isometric copies of some prototiles T^. 

Example: A diagram of the tile-substitution for the Penrose tiling [11] is shown in Figure 
[TJ Here, we use triangles as prototiles rather than rhombi or the celebrated kites and darts. 
Since the two prototiles are mirror-symmetric, but their substitutions are not, we have to 
respect the orientation of each tile. This is achieved by a black dot on each triangle. As 
indicated in the figure, iterating the substitution several times fills larger and larger portions 
of space, eventually resulting in a tiling of the plane. However, the proper definition of a 
tile-substitution and a substitution tiling is as follows. 

Definition 1.1. A (self-affine) tile-substitution is defined via a set of prototiles and a map a. 
Let T l: T 2 , ■ ■ ■ T m be nonempty compact sets — the prototiles — in U. d , such that cl(int(T )) = 
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Tj for each i < m. Let Q be an expanding linear map such that 

m 

(1) QT j = [jT i + V ij (l<j<m), 

i=i 

where the union is not overlapping (i.e., the interiors of the tiles in the union are pairwise 
disjoint); and each T>^ is a finite (possibly empty) subset ofW 1 , called digit set. Then 

<r(Tj:= {T i + V ij \i = l...m} 

is called a (tile-)substitution. By cr(Tj + x) := c(T-) + Qx and <r({T, T'}) := {a(T),a(T')}, 
a extends in a natural way to all sets {T^ + x k \ k £ /}. 

The sum T + = T + . . . , d ijfe } means {T + d^ x , . . . , T + dy fc }, with the convention 
T + = 0, if is empty. The set / in the last line of the definition can be any index set. 
In particular, if / = N, such a set can describe a tiling, given that the tiles do not overlap, 
and the union of all tiles is the entire space Before stating the definition of a substitution 
tiling, let us make two remarks. 

Remark: Often, a tile-substitution is defined without the requirement (pQ). To be precise, if 
the support of cr(Tj) equals QTj — that is, if ([I]) holds — , the substitution is called self-affine. 
Moreover, if Q is a homothety (i.e., just a multiplication by a scalar factor A), the substitution 
is called self-similar. In this case, A is called the substitution factor of the substitution. E.g., 
the expanding linear map for the Penrose tiling is Q = (jj), where r = (y/E + l)/2 is 
the golden mean. Thus, the substitution factor of the Penrose substitution in Figure [T] is r. 
Whenever the term 'self-similar tiling' ('self-affine tiling') is used in the sequel this is to be 
understood as 'self-similar substitution tiling' ('self-affine substitution tiling'). 

Remark: Every symbolic substitution (see Section ll.ip gives rise to a self-similar tile- 
substitution in R 1 in a canonical way, just by replacing each letter by an interval. This 
can be done in such a way that the tile-substitution is self-similar. Namely, one uses the 
entries of a positive PF-eigenvector of the substitution matrix, see [16], [8]. 

Definition 1.2. Let a be a tile-substitution. A tiling T is called a substitution tiling (with 
substitution a), if each finite subset of T is a translate of some subset of some a k {T-), where 
k > 0, and T- some prototile. The family of all substitution tilings with substitution a is the 
tiling space X CT . 

This definition of a tiling space is equivalent to the usual one [22], whenever a is primitive. 
Primitive means that there is k > such that a power of the substitution matrix (or 'incidence 
matrix') M a = (\T' i j\) 1<i j <m is strictly positive. Here, l^-l denotes the cardinality of D^. It 
is also consistent with the usual definition of the associated dynamical system ( X a , S) arising 
from a, see [16], as well as with our definition of X a in Section [1 .11 

1.3. Iterated function systems. The prototiles of a self-affine substitution can be derived 
from Q and (T> i j) 1<i - <m uniquely as follows: Multiplying the entire equation system ([TJ by 
Q^ 1 yields the corresponding iterated function system (IFS) for the tiles of a: 

m 

(2) T 3 = [jQ- 1 (T l +D !3 ) (l<j<m) 

i=l 
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Such an IFS is known to possess a unique nonempty compact solution (see for instance [5] 
and references therein). Since the prototiles of a fulfil ([2]), they are the unique solution of 
([2]). In particular, by Definition II .14 the IFS ([2]) is non-overlapping, that means, for each j, 
the interiors of the sets in the union are pairwise disjoint. Note that two parts Ti,Tj of the 
solution are allowed to overlap. The point is that ([2]) induces a partition of each Tj up to a 
set of measure 0. 

Originally, the term IFS was coined for the case m = 1 only. An IFS with m > 1 components 
is often called 'graph-iterated function system' or 'multi-component IFS'. Here, we use the 
term IFS for all of them. 

Remark: We just mentioned that each self-afnne substitution gives rise to a non-overlapping 
IFS ([2]). Conversely, by multiplying ([2]) with Q, each non-overlapping IFS with a common 
contracting map Q~ l defines a self-afnne substitution as in Definition 11.11 

Remark: An IFS as above can be written as a finite state automaton (FSA). The states are 
the T-. The input alphabet consists of all elements x in Ui<j j<m ^ij- The m P u t of a digit 
x changes the state from T- to T i; if x £ otherwise the state is changed to FAIL. For an 
example, see Figure El For details, cf. [23], [TO] . 

The following theorem ensures that it does not matter where the prototiles for a substitution 
a are located in space. The tiling space defined by a depends only on the way how tiles are 
dissected into smaller pieces, not where they are located. In particular, it does not matter 
whether the prototiles overlap or not. The important property is that the dissection of each 
prototile is non-overlapping. This theorem will be useful in the next sections, as well as for 
future work. 

Theorem 1.3. Let T> be the digit set matrix for an IFS in M d with a common contraction Q. 
If (A 1 ,A 2 , • • • , A m ) is a solution for V, then (A l + t 1 , A 2 + t 2 , ■ ■ ■ , A m + t m ) is a solution for 
{V l3 + (t. - QtA} x ^ m for all (t v t 2 , ...,tje (R d ) m . 

Proof. Throughout the proof, let (A±, . . . ,A m ) be a solution for V, i.e., 

m 
i=l 

First, let A'- = A - + 1 for some common t (1 < j < m). Then the approach 
A\ = QA[ +V u +x U • • • U QA' m + V ml + x 

A'm = QAi + T> lm + x U • • • U QA' m + T> mm + x 

yields 

A 1 +t = QA l + Qt + V n + x U • • • U QA m + Qt + V ml + x, 

and it follows x = t—Qt. It is straightforward to check that (A' l7 . . . , A' m ) = (A-^+t, . . . , A m +t) 
is a solution for {T>ij + (t - Qt)}i<jj< m - 
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Next, let t x = t 2 = ■ ■ ■ = t k = 0, t k+1 = ■ ■ ■ = t m = t and A'j = A- + t- for k + 1 < j < m. 
Then the approach 

A i = UtiQA + ^i u U? =J+1 QA' t + v hl + x 



(3) 



A' k+1 = UtiQA + ^k+x + y u l)Z k+ iQA + 'D i , k+1 + - 



yields y = t,x = —Qt and z = i — Qi. With these values, it is straightforward to check that 
(A x , . . . , A k , A k+1 + t, . . . A m + t) is a solution of ([3]). 

By combining these two results, we obtain the digit matrix for (A 1 +t 1 , A 2 + t 2 , ■ ■ ■ , A m + t m ) 
successively as follows. 

p(l) =X)+(t 1 -Qt 1 ) = Pii+Cii-Qti)}^^ has the solution (A 1 +t 1 ,A 2 +t 1 ,...,A m +t 1 ). 



T>W = V + (t 1 -Qt 1 ) + 



/ tj-tj ■■■ ^-^ 

-QCta-tj (B-Q)(* 2 -*i) ••• (-E-Q)(* 2 -*i) 



has the solution + A 2 + t x + (t 2 — ii), . . . , A m + t x + (t 2 — ii)) (where we use © with 
y = t = t 2 — t 1 , and 12 denotes the identity matrix), and so on, until 

/ tm-* m -l 



p(m) _ p(m-l) _j_ 



has the solution (A x + t l5 A 2 + t 2 , . . . , An-i + i m -l> An + *m-l + (* m - * m _i))- 

Now, 2?( m ) is of the desired form, which proves the claim. □ 

Regarding the action of a common regular linear map T on the solution (A\, . . . , A m ) with 
respect to T> rather than the action of different translations is much simpler. The same is 
true if all translation vectors t{ in the previous theorem had been equal. 

Lemma 1.4. Let T> be the digit set matrix for an IFS in M. d with a common contraction 
Q. Let T : R d -» R d be any invertible linear map. If (A l} A 2 , ■ ■ ■ , A m ) is a solution for T>, 
then (TA 1 ,TA 2 , . . -,TAJ is a solution for T(V), where T{V) = (T^)) ^..^ T{V t] ) = 

{ Td ijk\ d ijk £ v iji- 
Proof. 

m mm 

Aj = \jAi + Vij => TAj =T^\jAi + D y ) = \jTA l + T(V tj ). 

i=l i=l i=l 



□ 
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1.4. Model sets. In the sequel we consider self-similar tilings which arise from model sets. 
Such tilings are known as cut-and-project tilings [I], [H], [15]; see also [16], where the cut- 
and-project method for one-dimensional tilings is called 'geometric realization'. 

Definition 1.5. A model set is defined via a cut and project scheme, i.e., a collection of 
spaces and mappings as follows. 

Let G, H be locally compact Abelian groups, A be a lattice in G x H (that is, A is a cocompact 
discrete subgroup of G x H ), 7Ti : G x H — > G, 112 : G x H — > H be projections, such that 7Ti|a 
is injective, and ^(A) is dense in H. Let W C H be a compact set — the window — such 
that the closure o/int(W) equals W . 



G 




H 


U 


u 


U 


V 


A 


w 



Then 

V := {tti(x) \x € A, tt 2 (x) £ W} 

is called a model set. 

If fj,(d(W)) = 0, then V is called regular model set. 

The star map is the map * : 7i"i(A) — > H, x* = ^(^i 1 {%)) ■ 

We should mention that in many cases in the literature holds G = R d and H = R e . One may 
assume this reading this article without loosing much information. However, the notion of 
star-duality is tailored to the general case, and this is where we leave the concept of Galois- 
duality. There are cut-and-project tilings (in fact, a lot of them) of the line or the plane, 
where H is the field Q p of p-adic integers, or a product of those with some Euclidean space 
R e , see [3], [H], [H], [20]. 

In the context of symbolic substitutions [16], the setting is very much the same, but unfortu- 
nately the terminology differs. For instance, the cut-and-project scheme is called geometric 
realization, the window is usually called (generalized) Rauzy fractal etc. 

A model set is, by definition, a point set. There are several ways to assign a tiling to it. For 
instance, one may take the Voronoi cells of the point set as tiles. Whenever we can obtain a 
tiling T from a model set in such a way we call T a cut-and-project tiling. To be precise: If 
it is possible to assign a tiling T to a model set V such that V and T are mutually locally 
derivable (MLD) in the sense of @], we call T a cut-and-project tiling. Roughly spoken, MLD 
means that one can obtain V from T by local replacement rules and vice versa. 

The self-similar cut-and-project tilings are the objects we will consider in the following. We 
should mention that there is a huge number of such tilings, including the Fibonacci tilings 
and the Tribonacci tilings [16], the Penrose tilings and the Ammann-Beenker tilings [IT], the 
Chair tilings and the Sphinx tilings [3] . A collection of several additional examples is available 
online [9]. 

2. The *-dual substitution 

It is known that the window of a self-similar model set is the solution of an IFS, see for 
instance [12], [H], [21] ■ Since each non-overlapping IFS defines a tile-substitution (see the 
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remark after ([2])), several authors considered this 'dual substitution tiling' (or better, the 
substitution) induced by the IFS of the window. In [23], Thurston uses a different construction 
based on Galois conjugation in the field Q(A), where A denotes the substitution factor of a 
self-similar tiling. One can utilise Theorem 11.31 to show that both approaches yield the same 
tiling spaces. Thurston called the tilings arising from this construction 'Galois duals' of the 
original tilings. We proceed by stating the definition of the *-dual tiling of a self-similar 
cut-and-project tiling, which is a generalization of Galois-duality. 

In general, the groups G and H in Definition 11.51 can be arbitrary locally compact Abelian 
groups. However, dealing with substitutions and linear maps, we require G and H to be 
equipped with a vector space structure in the sequel. This is no rigid restriction, since all 
examples in the literature fulfil this condition. 

Furthermore, in this section we require the expanding linear map Q : G — > G to be as follows: 
Let T: GxH^GxHbea hyperbolic linear map (i.e., there is no eigenvalue e of T with 
|e| = 1) such that T\q is an expansion and T\h is a contraction. Let Q = T\q. Again, this is 
not really a restriction, the majority of well-studied tilings fit into this framework. 

Definition 2.1. Let T be a self-similar cut-and-project tiling with substitution a and expan- 
sion Q. Let V = (Pij)i<i,j<m be a corresponding digit matrix as in Definition \1.1\ such that 
a is the substitution induced by V. Furthermore, let Qtti(A) C 7Ti(A), and each d G V- be 
in 7q(A), with tt\ and A as in Definition 11.51 Then (T>*) T defines a substitution a* with 
expansion Q' := (T\h)~ 1 , which is called the *-dual substitution. 

Here, V* = (Z$)i<y < m > and V % = {d*\de V i3 }. 

If the substitution is one-dimensional and self-similar with substitution factor A, where A is a 
unimodular PF-number, then Q' = (diag(A 2 , A 3 , . . . , X N )) 1 , where A 2 , . . . A^ are the Galois 
conjugates of A, see [ID], [8]. Therefore the name 'Galois duals' was coined in [23] . 

Remark: Note that the star-dual substitution also defines self-similar cut-and-project tilings, 
with the roles of G and H in Definition 11.51 reversed. In particular, a star-dual tiling of some 
tiling in G is a tiling in H. There are tilings of the Euclidean plane (or the line) with p-adic 
internal spaces H, see [3], [20]. For instance, the star-dual tiling of a plane tiling may be a 
tiling in Q p x Q p , where Q p is the field of p-adic integers, which is a locally compact Abelian 
group as well. 

Remark: It follows immediately from the definition that (c*)* = a. By Lemma 11.41 we can 
freely choose between applying the star map to T> or The resulting substitutions are 

the same, up to a scaling of the prototiles. 

Example: (Fibonacci squared) Consider the word substitution C, on the alphabet {a, 6} given 
by £(a) = aab, ((b) = ab. The substitution matrix (or 'incidence matrix') is = \\\)- It 
is well known — and easy to see — that the PF-eigenvalue (the leading eigenvalue) of M is 
the substitution factor for the tile-substitution. Moreover, the entries of a left eigenvector of 
M, corresponding to the PF-eigenvalue, contains the relative lengths of the prototiles for a 
self-afnne tile-substitution [8], [16]. The leading eigenvalue in this case is r 2 = ^ 3 , which 
is the substitution factor of the desired substitution. Requiring the prototiles to be intervals, 
the lengths can be chosen as £ a = 1, £ b = t _1 . (Here is some freedom, we may choose any 
multiples of these lengths as well.) One appropriate tile-substitution is given by the digit set 
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Figure 2. The automaton for the IFS for the squared Fibonacci substitution 
and its dual automaton. 



matrix 

'{0,1} W 



{2} {!}, 

The prototiles T a = [0, 1] and T b = [0, r _1 ] arise as the unique compact nonempty solution of 
the corresponding IFS 

T a = r- 2 T a U r-\T a + 1) U r-\T b + 2), 
T b = T' 2 T a UT-\T b + l). 



We proceed by considering the digit set of this IFS, V = ( ^?,' r _2i^ r^-n )> m agreement with 



{2r-2}{^ 2 }. 

the last remark, and in order to illustrate the action of Galois conjugation. The IFS (j3]) can 
be represented by the automaton in Figure [21 left. (Vice versa, the automaton gives rise to 
the IFS ([!])). The dual substitution arises from the IFS defined by (T>*) T , or equivalently: 
from the dual automaton, where the orientation of each edge is reversed and the star map is 
applied to each edge label, see Figure El right. In this case, the star map is uniquely defined 
by 1* = 1, t* = —t~ 1 (each r which occurs is replaced by its Galois conjugate — r _1 ). Thus, 

( r - 2 )* = (1 - t- 1 )" = (2 - r)* = 2 + t- 1 = 1 + r = r 2 . 

The corresponding digit matrix, defining the dual substitution £*, is 



*nT_/{0,t 2 } {2r 2 } 



{0} {r 2 } 



The solution of the corresponding IFS — with digits r 2 d it - k , where d^ k are the digits in 
(V*) T — yields the prototiles [0, r] and [r, r + 1] for £*. Since this dual tile substitution is 
nice in the sense that the tiles are intervals — rather than some disconnected, fractally shaped 
sets — it can be written as a word substitution a* t— > a*b*a* , b* ^> b*a* . In particular, even 
though C 7^ C^i they generate — up to a scaling by r — the same tiling spaces (the same 
tilings), resp. the same biinfinite words. 



3. Dimension 2 

In dimension one there is a canonical way to associate a tiling to a model set, which is, 
by definition, not a tiling but a uniformly dense point set: The points define a partition of 
the line into intervals, these intervals are the tiles. Analogously, a tiling of M 1 by intervals 
defines a point set by considering the set of all vertices, i.e., the set of all boundary points of 
the intervals. In particular, by assigning each interval to its left boundary point, there is a 
one-to-one correspondence between tiles of the tiling and points in the point set. 
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In two dimensions the situation becomes more difficult, at least if one wants to study the 
different types of tiles/points x through their position x* in the window W. It is not obvious 
which way to identify tilings with point sets (and vice versa) is the best one, or the most 
natural one. 

On one hand, it seems natural to consider the vertices of a tiling T. In fact, this point of 
view is taken in many cases, see e.g. [2], [12]. Then, in general, there is no natural one-to-one 
correspondence between tiles and points. In particular, the positions x* in W, where x is 
some vertex, may give us no relevant information about the tiles in T. One simple reason is 
that, in general, the number of tiles per unit area differs from the number of vertices per unit 
area. 

On the other hand, in order to assign a point to each tile, one can use control points |22| . 
Then there is a one-to-one correspondence between tiles and points. But there are many 
ways to choose such control points, thus one may end up with different windows for different 
choices. 

In the sequel we avoid all the problems mentioned above by computing some *-dual substi- 
tutions of prominent two-dimensional tilings. These computations follow the methods in [10] 
closely. In the framework of Definition II. 1[ the number of prototiles of a two-dimensional 
substitution is considerably large, since we identify prototiles only up to translation, not up 
to isometries. For instance, the Penrose tiling uses only two different triangles (up to isome- 
tries) as prototiles, see Figured! But each of those occurs in 10 different orientations in the 
tilings. Moreover, we need to distinct a tile and its mirror image (see the dots in Figured]), 
thus there are 40 prototiles altogether with respect to translations. Rather than working with 
a 40 x 40 digit matrix, we will slightly change our framework: Instead of digits representing 
translations, we go over to 'digits' which are isometries. The entries in the according matrices 
now are functions rather than translation vectors. Each function is of the form x QRx + d, 
where Q is — as usual — an expanding linear map which is the same for all functions, R 
is some rotation or reflection, and d is a translation vector, as above. In other words: We 
consider matrix function systems instead of digit matrices, cf. [13]. Since we are considering 
plane tilings, we will identify M 2 with the complex plane C. Then, a rotation about the origin 
is just a multiplication with some complex number z with \z\ = 1. Each reflection in the 
complex plane can be expressed as complex conjugation, followed by some rotation. The 
expanding map Q will always be a homothety in the remainder of this section. Our goal is to 
express all translations and rotations, and thus all maps f i occurring in the IFS under con- 
sideration, in the ring Z[£], where £ is some root of unity. To obtain the dual substitution, we 
consider the inverse maps ff 1 and apply a primitive element of the Galois group GaZ (Z[£]/Z) 
to them. In plain words: we replace each £ occurring by an appropriate Galois conjugate 
£' in order to obtain the desired maps (ff 1 )*- This is in fact how the star map acts in the 
following cases, compare [2]. To simplify notation, we denote (Z -1 )* by fK 

3.1. The *-dual of the Penrose tiling. Consider the prototiles of the Penrose tiling as in 
Figure [Tl It is known that the vertices of a Penrose tiling are closely related to the cyclotomic 
field Q(£ 5 ), where £ := £ 5 = e~ . E.g., all vertices, and all edge lengths, can be expressed 
in the ring Z[£] of integers in Q(£). Thus, let the small triangle S have vertices 0, —1 and 

27ri3 a liri q 

e io = — £ , and the large triangle have vertices 0, r and eTo~ = — £ d . Then, the IFS for these 
tiles — arising from the substitution — uses the common contracting factor r _1 = £ + £ 4 , 
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Figure 3. The automaton of the IFS of the Penrose tiles (left), compare ([5]), 
and its dual automaton (right). 



and the IFS reads: 

(5) S = f x {S) U f 2 (L), L = f 3 (L) U / 4 (fi) U f 5 (L), 

where 

fi(x) = (C + e)(-^)x-l-f, 

f2(x) = (t+z 4 )(-ox-i-e-z\ 

/ 3 ^) = (e+c 4 )e 3 x + i+c, 

h(x) = (t+e)(-e)x+i-e, 

h(x) = (C + e)(-l)x + l-f-f. 

The corresponding automaton is shown in Figure [3] (left). Note, that in [10] the maps are 
different and do not define the Penrose substitution. Thus, the author of [10J did not compute 
the *-dual (or 'Galois-dual') of the Penrose substitution, but of a different one, using the same 
triangles S and L as prototiles. 

In order to calculate the *-dual substitution, consider the dual automaton (see Figure [3j 
right), resp. the dual IFS. Let <p G Gal(Z[£\) with ip(£) = £ 3 . Note that <p is uniquely defined 
by this requirement. The /? are obtained by applying 92 to ff : 



/!(*) 


= (Z + Z 4 )(fx 


+ f+f) 


fkx) 


= (t + Z 4 )(fx 


+ c + e 2 + e 4 


fkx) 


= (c + e 4 )(-^ + e+^ 4 ) 


fkx) 


= (£ + e 4 x^- 


-t+t 2 ) 


fkx) 


= (c + e 4 )(^- 





The corresponding dual IFS reads: 

(6) S* = /{(fif*) U /1(L*), L* = /»(£*) U /J(L*) U /«(L*). 

A numerical computation of the solution yields the image in Figure 0] (left). In the figure, 
the two components of the solution are torn apart; actually, the solutions do overlap. This 
fact is not relevant, due to Theorem 11.31 The important point is that the dissection of each 
prototile is non-overlapping. 

It is reasonable to assume that the solution consists of two triangles similar to the Penrose 
triangles. In order to make this precise we list the coordinates of the triangles: 

l* : i + e 2 + e 4 , i+c + e 3 , i+f + f, s* : ^ 2 ,e 3 ,2e 2 -e + e 4 -2e 3 . 
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Figure 4. A numerical computation of the solution of the dual IFS of the Pen- 
rose substitution (left), and the Tubingen triangle substitution (right). One 
may guess from the images that this is the *-dual of the Penrose substitution, 
which turns out to be true. 




Figure 5. The substitution rule for the Ammann-Beenker tiling (left), with 
prototiles S, L, and the star-dual substitution (right), with prototiles S*,L*. 
Even though the two substitutions look different on a first glance, they define 
very similar tilings, see Figured! 

All this information provided, it is an easy — but pretty lengthy — exercise to check that these 
triangles are in fact a solution of ([6]). The tile-substitution arising from ([6]) is well-known: It 
is the Tubingen Triangle tile-substitution (see Figured! right). See [9] for more details and 
images of this substitution, as well as several others. Altogether we have established: 

Theorem 3.1. The -k-dual of the Penrose substitution is the Tubingen triangle substitution, 
and vice versa. 

3.2. The *-dual of the Ammann-Beenker tiling. In a very similar fashion, one can 
compute the star-dual of the Ammann-Beenker tiling [11]. The substitution rule of the 
Ammann-Beenker tiling is shown in Figure [5] (left). A part of an Ammann-Beenker tiling is 
shown in Figured] (left). The Ammann-Beenker tiling is closely related to the cyclotomic field 
Q(£ g ), where £ := £ 8 = e~ . Its *-dual can be calculated in a very similar fashion as in the 
case of the Penrose tiling. The substitution factor is A = 1 + y/2 = l+ £ + £ 7 = l + £ + £ _1 , 
its inverse A -1 = —1 + £ + Let S be the triangular prototile, with vertices 0, i, 1 + i, and 
let L the rhombic prototile with vertices 0, 1, 1 + The IFS for these prototiles is 
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where 





MaC *>k 

k <Jk ^kHUk 

rjarir vhif 
iS^BlS^^&Jk ^Bk v&SMifa 



.^k .^k fiSSffi? 



kHUk aQa 
k ^Ok *w ^k 



Figure 6. Patches of the Ammann-Beenker tiling (left) and its star-dual 
(right). The patches are equivalent in the sense, that one obtains the right 
patch by dividing each tile in the left patch into four tiles. 




L = 



~ fi(L) U f 2 (L) U / 3 (L) U f & (S) U f 7 (S) U f 8 (S) U f 9 (S), 
S = ML) U ML) U f w (S) U / U (S) U f 12 (S), 



Mx) = \- i ex + l + t 

Mx) = \-^x+i+^+e+e, 

'-^ 7 x + l + £ + f, 

'x+z+e, 



Mx) = A _1 a 

/ 5 (x) = A- 1 e 2 x + ^, 

/ 6 (x) = A- 1 e 2 x + ^, 



/ 7 (x)=A- 1 C 3 x + l + 2C + C 2 + f, 

/ 8 (x) = A- 1 c 6 x + 2 + e+^ 2 + e 7 , 

/ 9 (x) = A-^ 7 x + 1, 

Mx) = \- i fx+z+e, 
f u (x)=\- i fx+i+z+e, 
f l2 ( x ) = \- i f x+ 2t+e. 



The functions ff of the dual IFS are easily obtained by considering the inverse functions /~ , 
and replacing each £ occurring by its Galois conjugate £ 3 (or £ 5 or £ 7 , the resulting dual tiling 
is always the same). The dual IFS reads: 

V = fl(L*) U /{(£*) U /*(£*) U fl(S*) U 

S* = fl(L*) U /J(L*) U fl(L*) U /*(£*) U /f (S*) U fl(S*) U /*,($*). 

Its solution are the two tiles shown in Figure [5] (right): S* is an isosceles orthogonal triangle, 
and L* is a thin orthogonal triangle. Note that the original triangular prototile S can be 
dissected into two copies of S*, and the rhombic prototile L can be dissected into four copies 
of L*. In fact, unlike in the Penrose case, the *-dual tiling of the Ammann-Beenker tiling is 
closely related to the Ammann-Beenker tiling itself, compare Figure [6) dissecting each L into 
four tiles, and each S into two tiles, yields the *-dual tiling of the Ammann-Beenker tiling. 
Altogether we have established: 
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Theorem 3.2. The -k-dual of the Ammann-Beenker tiling is mutually locally derivable with 
the Ammann-Beenker tiling. 

4. Self-duality 

The last example, the Ammann-Beenker tiling and its *-dual, motivates the question whether 
there are substitutions which are strictly self-dual. 

Definition 4.1. Let a be a tile- substitution such that X CT consists of cut-and-project tilings. 
If X CT = X CT * , then a is called self-dual (with respect to ^-duality) . 

Thus, the example in Section [U Fibonacci squared, is self-dual. Since the star-dual tiling 
of a tiling in G is a tiling in H (cf. Definition II. 5p . a trivial necessary condition for a tiling 
to be self-dual is G = H. Hence, if the substitution factor A of a self-dual substitution is a 
unimodular PV-number, then A must be a quadratic algebraic number. A refined necessary 
condition is obtained by considering the substitution matrix M a . 

Lemma 4.2. If a is a self-dual substitution in R , then the following property holds for the 
substitution matrix M a = 

(7) (M a ) T = PM a P~\ 

where P is some permutation matrix. 

Proof. The result follows immediately from M CT * = (|(Z> i;? -)*|) T , by definition of a*. The 
tiling spaces X<j defined by a substitution are unique up to enumeration of the tiles, thus the 
permutations enter. □ 

In fact, permutations has to be taken into account. There are self-dual substitutions where 
(M a ) T M a ; in other words: where M a is not symmetric. For instance, the dual substitution 
of o~ '. a — > ab, b — ► aab is a* : a* — > b*b*a*, b* — ► 6*a*. By identifying the first tile — a 
- of a with the second tile — b* — of the dual substitution a*, one obtains that the two 
substitutions define the same tiling spaces. Hence a is self-dual. The substitution matrix 
M a = (21) is n °t symmetric, but 

thus © holds. 

A sufficient condition for a substitution to be self-dual is an immediate consequence of The- 
orem [L3j 

Corollary 4.3. Let a be a substitution with digit matrix T> and a* its dual substitution, with 
digit matrix V = (T>*) T . If there are vectors t 1 ,t 2 , ... ,t m € M. d , A := (t - — Qtj)\<ij< m and 
a permutation matrix P such that 

PV'P- 1 =T> + A, 

then a is self-dual. □ 
This result gives a purely algebraic test for self-duality. 
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5. Outlook 

The notion of star-duality is a concept tailored to arbitrary model sets in any locally compact 
Abelian group. Restricted to the unimodular case, it coincides with other notions of duality, 
in particular the 'natural decomposition method' in |21j . and the dual tilings considered by 
P. Arnoux, S. Ito and others, see for instance [7], [T7|. This is the reason why we were 
lazy at some occasions, using the term 'self-dual' instead of 'self-*-duaF or 'self-dual w.r.t. 
star-duality'. Theorem 11.31 can be utilised to show the equivalence of the concepts mentioned. 

The benefit of star-duality seems to be its generality, and its accessibility by purely algebraic 
methods. This paper tries to make first steps in this direction. For future work, it should be 
promising to apply star-duality to the case of word substitutions on two or three letters. In the 
case of Sturmian substitutions (symbolic substitutions on two letters, where the corresponding 
biinfinite words contain exactly n + 1 subwords of length n, for all n > 0, see [16]) the results 
in [6] can be used to give a complete characterisation of self-dual (Sturmian) substitutions. 
For details we refer to future work. 

A further question is: are there self-dual substitutions in dimension d > 2? The Ammann- 
Beenker tiling is very close to be self-dual, but is is not exactly self-dual. In dimension 
one, there are several self-dual substitutions. Trivially, these examples generalise to higher 
dimensions, by using the Cartesian product (for instance, Fibonacci times Fibonacci). But it 
remains to find a two-dimensional self-dual substitution apart from that. 
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